Incommensurate spin density wave in Co-doped BaFe 2 As 2 
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57 Fe Mossbauer spectroscopy measurements are presented in the underdoped Ba(Fei_a;Coa;)2As2 
series for 1=0.014 (T c < 1.4 K) and a;=0.03 and 0.045 (T c ~ 2 and 12 K respectively). The spectral 
shapes in the so-called spin-density wave (SDW) phase are interpreted in terms of incommensurate 
modulation of the magnetic structure, and allow the shape of the modulation to be determined. In 
undoped BaFe2As2, the magnetic structure is commensurate, and we find that incommensurability 
is present at the lowest doping level (a;=0.014). As Co doping increases, the low temperature 
modulation progressively loses its "squaredness" and tends to a sine-wave. The same trend occurs 
for a given doping level, as temperature increases. We find that a magnetic hyperfine component 
persists far above the SDW transition, its intensity being progressively tranferred to a paramagnetic 
component on heating. 

PACS numbers: 74.70.Xa, 75.30.Fv, 76.80.+y 



I. INTRODUCTION 



Investigation of the magnetic properties of the newl 
discovered layered Fe-based pnictide supercondutors [l 
with a microscopic probe such as Mossbauer spectroscopy 
is of interest for (at least) two reasons: i) to determine the 
characteristics of the expected spin-density wave (SDW) 
order and ii) to give experimental evidence concerning 
the problem of the interaction and/or local coexistence 
between superconductivity and magnetic ordering of the 
Fe moments. In the so-called 122 family with Ba, the par- 
ent (non superconducting) compound is BaFe2As2 and 
doping with Co substituted for Fe has been shown to in- 
duce superconductivity at low doping levels @. In the 
Ba(Fei_ x Co x )2As2 series, the phase diagram as a func- 
tion of Co concentration x is "bell-shaped", with T c (x) 
presenting a maximum of 26 K for £=0.07 0, H[. Be- 
low this optimal concentration, the material shows ap- 
parently the persistence of an antiferromagnetic (AF) 
structure in the superconducting phase, the AF or spin 
density wave (SDW) transition temperature decreasing 
rapidly as x increases. Neutron diffraction in the un- 
doped compound BaFe2As2 [I| has determined the AF 
magnetic structure below TV ~143K to be a priori com- 
mensurate. 57 Fe Mossbauer spectroscopy [f| 0] and 75 As 
NMR [8] in BaFe2As2 have confirmed the commensura- 
bility of the magnetic structure. 

In this work, we focus on the first problem evoked 
above: we performed a 57 Fe Mossbauer study of 
Ba(Fei_ x Co x )2As2 in the underdoped region in order to 
investigate the SDW state. In our samples with Co dop- 
ing £=0.014, 0.03 and 0.045, the resistivity data show an 
inflexion point [H respectively around 115 K, 92 K and 
64 K (see FigJTJ). According to neutron diffraction data 
in materials with similar Co dopings 0-[H|, this tem- 
perature corresponds to the SDW-paramagnetic transi- 
tion. We find that the spectra in the magnetic phase 
have an unusual shape, which we interpret as due to 
an incommensurate modulation of the AF structure, i.e. 
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FIG. 1: [color online] Resistivity vs. temperature for 
BaFe 2 As 2 with Co dopings x=0.014, 0.03 and 0.045 from 
Ref.0]. The data for £=0.014 have been shifted by 100 pfl.cm 
for clarity. The full arrows mark the inflexion point of p(T), 
attributed to the onset of magnetic ordering (Tsdw), and the 
dash-dotted arrows the superconducting transition tempera- 
ture. 



to an incommensurate (IC) spin density wave. We also 
find that the magnetic hyperfine subspectrum persists far 
above Tsdw, i- e - that AF correlations are present up to 
Taf ~ 1-5 Tsdw- 



II. EXPERIMENTAL AND SAMPLE 
CHARACTERISATION 

The samples are single crystals prepared by the self- 
flux method, as described in Ref.[4|. It was checked 
that Co doping is uniform and that superconductivity 
occurs in the bulk [13, LL3J . For the Mossbauer measure- 
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TABLE I: In Ba(Fei_ : 
and maximum value h r , 
tion at 4.2 K . 



Coa;)2As2: Fourier coefficients /i2fc+i 
i (in T) of the hyperfine field modula- 
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ments, the crystals were finely ground. Powder samples 
made by solid state reaction were also prepared, with 
x=0 and 0.03, for comparison purposes. The spectra 
were recorded using a commercial Co*:Rh 7-ray source 
mounted on an electromagnetic drive with linear velocity 
signal. Room temperature spectra ensured that no trace 
of FeAs (with a quadrupole splitting AZ?Q=0.69mm/s) 
was present. 

In the high temperature tetragonal phase (space group 
IA/mmm), the Ad Fe site has fourfold symmetry (4m2). 
At low temperature, slightly above Tsdwj a small or- 
thorhombic distortion takes place (space group Fmmm), 
and the local symmetry at the 8/ Fe site is lowered (222). 



III. THE 57 FE MOSSBAUER SPECTRA AT 4.2 K 

The spectra at 4.2 K in the three investigated com- 
pounds and in undoped BaFe2As2, are represented in 
Fig. IIII1 In BaFe2As2, one observes a single magnetic 
hyperfine pattern, in agreement with Ref.j?;], with a hy- 
perfine field H h f = 5.4(1) T. By contrast, the spectral 
shapes in the Co-doped compounds are rather far from 
this standard 6-line pattern. Saturation effects can be 
discarded since the resonant absorption is modest; for 
x=0.03, where it amounts to 7%, we checked that an ab- 
sorber with 2% absorption yields exactly the same spec- 
trum. Texture effects are also negligible since a powder 
absorber (for x=0.03) yields spectra quasi-identical with 
those of the single crystal sample. The high statistics 
reached in these experiments allows us to distinguish fine 
details of the spectral shape, which is a prerequisite for 
the fits to be described below. 

We interpret these unusual spectral shapes as due to a 
distribution of hyperfine fields arising from an IC mod- 
ulated magnetic structure. Indeed, in this case, the 
57 Fe nuclei are submitted to hyperfine field values rang- 
ing from zero to a maximum value, and although the 
Mossbauer data cannot give access to the propagation 
vector q, the shape of the distribution reflects the par- 
ticular shape of the modulation. In order to describe 
the influence of the moment modulation on the spectral 
shape, we assume a collinear magnetic structure, with 
no ferromagnetic component, like in undoped BaFe2As2 
where the Fe moments lie along the orthorhombic a axis 
p|. We also assume the hyperfine field is collinear with 
the Fe moment (see discussion below). Then, along the 
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FIG. 2: Mossbauer absorption spectra on 57 Fe at 4.2 K in 
Ba(Fei_ a; Co ;c )2As2 for x=0, 0.014, 0.03 and 0.045. Except for 
x=0, the solid lines, which are masked by the data due to the 
goodness of the fits, are simulations with an incommensurate 
modulation of hyperfine fields (see text). 
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FIG. 3: [color online] Modulation of the hyperfine field along 
the direction of q at 4.2 K in Ba(Fei_ a; Co a ;)2As2 for £=0.014, 
0.03 and 0.045, derived from the Mossbauer spectra. Insert: 
variation with the Co concentration x of the maximum hyper- 
fine field value (circles); for x=0, hyperfine field in BaFe2As2 
(square). The dashed line is a guide for the eye. 
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direction OX parallel to the propagation vector q, the 
modulation of the hyperfine field is defined in terms of a 
Fourier series as: 



H hf (qX)= h 2k+i sin[(2fc + l) 5 X], 



(1) 



where the hi are the odd Fourier coefficients of the mod- 
ulation. The spectra are then fitted to a superposition 
of individual hyperfine 6-line patterns with discrete f 
values according to Eq. (JXJ) with a linear mesh along a pe- 
riod (0< qX < 27r). The number of Fourier components 
in the fit is increased until a good reproduction of the 
lineshape is obtained. At 4.2 K, the very good fits shown 
in Fig lllll were performed with 6, 5 and 3 Fourier com- 
ponents for x=0.014, 0.03 and 0.045 respectively. The 
resulting hyperfine field modulations are shown in Fig|3l 
and the values of the Fourier coefficents are given in Ta- 
ble HI The isomer shift at 4.2 K with respect to a-Fe is 
0.54(1) mm/s and a very small quadrupolar interaction 
with AEq ~ —0.04(2) mm/s is needed to reproduce the 
slight asymmetry of the spectra. 

The 4.2 K modulation is seen to depart more and 
more from "squaredness" (which occurs for x=0) as x 
increases; for cc=0.045, the modulation is close to a pure 
sine-wave. This is reflected in the number of Fourier co- 
efficients needed to reproduce the lineshape (see TablelJ), 
which decreases as x increases. The maximum value of 
the hyperfine field at 4.2 K decreases as x increases (see 
insert in Fig©. For £=0.045, it is 75% of the £=0.014 
value. 



IV. THERMAL VARIATION OF THE 
MOSSBAUER SPECTRA 
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FIG. 4: [color on line] Mossbauer absorption spectra on 57 Fe 
at selected temperatures in BaFe2As2 with x=0.014 Co dop- 
ing. The lines are fits to an incommensurate modulation 
of hyperfine fields; the red subspectrum at 130 K is a single 
Lorentzian-shaped line. 



The thermal variation of the spectra for £=0.014, for 
which T S dw —115 K (FigO}, is shown in FigED No im- 
portant change occurs up to about 80K, then the mod- 
ulation changes shape, tending towards a pure sine-wave 
as temperature is further increased. A striking feature 
is observed when comparing Fig lllll and Fig llVI increas- 
ing the Co doping has the same effect on the modulation 
at 4.2 K as increasing the temperature for the £=0.014 
sample. The maximum value of the hyperfine field shows 
little thermal variation, as well as the main Fourier co- 
efficient hi. Above 115 K, a single line, characteristic of 
the paramagnetic phase, is present in the spectra and it 
grows on heating at the expense of the magnetic sub- 
spectrum (see red line in the 130 K spectrum in Fig UVjl . 
The latter actually persists in a rather large temperature 
range above Tsdw> up to 150 K at least. In this temper- 
ature range, very good fits of the spectra were obtained 
with two components: a single Lorentzian-shaped line 
and a magnetic subspectrum with a modulation identical 
to that of the 115 K spectrum. For this latter component, 
a common scaling factor of the Fourier coefficients was 
introduced in order to allow for a decrease of the modu- 
lation amplitude on heating; it remains above 0.9 up to 
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FIG. 5: [color online] In BaFe2As2 with x—O.OlA Co doping, 
thermal variation of the fraction of the single line component 
in the Mossbauer spectra, for increasing and decreasing tem- 
perature runs (left scale), and of the resistivity (right scale). 
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FIG. 6: [color online] In BaFe2As2 with x=0.03 Co doping: 
Top thermal variation of the main Fourier component hi of 
the hyperfine field modulation (the line is a guide for the eye) 
and of the 4 next components /13 — /19 (insert) below Tsdw = 
92 K; Bottom thermal variation (on heating) of the fraction 
of the single line (left scale) and of the resistivity (right scale). 



not be distinguished, but the coexistence region probably 
extends up to 130-140 K. 

For 2=0.045, the modulation remains unchanged up 
to about 60 K, close to Tsdw, then a single line grows in 
the spectrum on further heating. However, due to the 
lack of resolution of the spectra, no quantitative assess- 
ment of the fraction of single line can be done. For this 
doping level, the spectrum at 4.2 K belongs to the super- 
conducting phase, but we postpone a discussion of the 
mixed phase (superconducting/SDW) to a future publi- 
cation. 

According to the neutron diffraction data [3-[Tl| ■ the 
transition to a long range ordered SDW state occurs at 
Tsdw- Actually, this transition is not "seen" in the spec- 
tra of our underdoped samples; rather, above Tsdw, a 
paramagnetic component (single line) appears and co- 
exists with the magnetic hyperfine component. The 
latter progressively dwindles as temperature increases, 
its (maximum) hyperfine field decreases slowly on heat- 
ing, and so does the Fe moment (see discussion below). 
We also observe a small hysteresis (for x=0.014). This 
points to the persistence of AF correlations, probably 
short range and dynamic, far above Tsdw- We estimate 
that a "fully paramagnetic" phase is recovered around 
Taf ~ 1.5 Tsdw- This behaviour could be linked with 
the strong bidimensional character of the magnetic inter- 
actions in the Fe layers [HI, [l6[ . 

As to the structural transition which, in the Co doped 
compounds, occurs at a temperature higher than Tsdw 
S El) nothing can be said from the Mossbauer data 
due to the very weak quadrupolar interaction in these 
compounds. 



the highest temperature. The thermal variation of the 
fraction of the single line derived from these fits is shown 
in Figj5] a small hysteresis is observed and a unique line 
is recovered only above 150-160 K. 

For the sample with £=0.03, the thermal variation of 
the 5 Fourier components was studied in more detail (see 
Fig|6] top) and a similar trend is observed. The main 
component hi is an order of magnitude larger than the 
higher harmonics, and its variation with temperature is 
weak: it falls by about 20% between 4.2 and 90 K, which 
is close to Tsdw (see Fig|TJ). The higher harmonics are 
constant up to 60 K, then they drop and /13 becomes neg- 
ative and almost opposite to /15 , the other two harmonics 
being smaller. In terms of the shape of the modulation, 
it does not change up to about 60 K, and on further heat- 
ing it evolves smoothly towards an almost pure sine- wave 
at 85 K. The spectrum at 85 K for x=0.03 is very similar 
to the 4.2 K spectrum for x=0.045, illustrating the above 
mentioned trend that increasing x at 4.2 K has the same 
effect on the shape of the modulation as increasing T for 
a given x. Above 90 K, i.e. above Tsdw, a coexistence 
region is present, like in the sample with x=0.014. The 
thermal variations of the fraction of the single line and 
of the resistivity are represented in FigjS] bottom. Above 
115 K, the two components of the Mossbauer spectra can- 



V. DISCUSSION 

A. Disorder induced by Co substitution 

Since these materials are (weakly) substituted com- 
pounds, one may wonder whether the observed line- 
shapes, where a distribution of hyperfine fields is clearly 
present, could be accounted for by some disorder induced 
spread in the hyperfine field values. In the layered struc- 
ture of the pnictides, a Fe atom has 4 in-plane nearest 
neighbors, and for Co dopings £=0.03 and 0.045, it re- 
sults that resp. 12% and 18% Fe atoms have one Co atom 
as nearest neighbour. Assuming that the hyperfine field 
is different for a Fe atom with 4 Fe and with 1 Co - 3 Fe 
as nearest neighbors, then the corresponding Mossbauer 
spectrum would show a two peak hyperfine field distri- 
bution with, for instance for x=0.03, relative weights of 
88 and 12%. The distribution associated with the 4.2 K 
modulation for x=0.03 is shown in FigJTJ it presents 3 
peaks, each with 25-30% relative intensity, which there- 
fore cannot be due to the above hypothesis of different 
environments of a Fe atom. The hyperfine fields of Fe 
atoms with 4 Fe or 1 Co - 3 Fe are probably very close 
due to the delocalised nature of the 3d electrons. 
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TABLE II: Saturated Fe moment m and hyperfine field Hhf, 
and the ratio r = H^i/m in AFe2As2 for A= Ba, Ca and Sr. 



2 3 4 

Hyperfine field (T) 



FIG. 7: [color online] Hyperfine field distribution associated 
with the hyperfine field modulation at 4.2 K for the £=0.03 
sample (see Table [I]). The positions of the peaks near 2.8 and 
4.3 T correspond to the quasi-plateaus of the modulation, and 
the (unphysical) divergence to the maximum hyperfine field 
where dHhf /d(qX)—Q (see Fig(3|. Insert: integrated weight. 



B. Other evidences for incommensurability 

An incommensurability of the AF magnetic order in 
Co-underdoped BaFe2As2 has been inferred by NMR 
measurements [l2l . [l4| from the anisotropy of the line 
broadening. The IC wave-vector was estimated in 
Ref.[l2| to depart only slightly (by a few percent) from 
the commensurate value q=(101) in the orthorhombic 
cell. /xSR data in a sample with x=0.04 [? ] can also be 
interpreted in terms of an IC magnetic structure. The 
neutron diffraction measurements [9l-lll| did not report 
any incommensurability, but the resolution in q-space is 
probably not sufficient for this purpose; we note that 
the (a,b) orthorhombic splitting could not be resolved in 
these experiments as well. In LaOFeAs, a calculation of 
the Lindhardt response function xo (<z) was performed for 
the undoped system and for the system with £=0.1 F 
doping [21(. It is shown that F doping shifts the maxi- 
mum of xo (q) from the M point of the Brillouin zone to 
an IC q vector, inducing thus an IC magnetic structure. 
The appearance of an IC SDW upon doping is also pre- 
dicted in Ref.[? ] In view of the similarities between the 
1111 and 122 families of pnictides, these results could also 
hold for the latter family, accounting for the fact that a 
very small doping level (x=0.014) is enough to push the 
material towards an IC SDW. 



C. The hyperfine constant and the Fe moment 

The saturated hyperfine field in BaFc2As2 is 5.4 T [7] 
and the maximum hyperfine field at 4.2 K in the Co 
doped compounds has comparable magnitude (see insert 
of FigJ3J). These are unusually small values for Fe, even 
in an itinerant magnet. In intermetallic compounds, it is 
usually assumed that the Fermi contact interaction with 
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s-electrons is the dominant hyperfine coupling, implying 
that the Fe moment m is proportional to the hyperfine 
field Hhf 22j. This is true in a-Fe, with a hyperfine con- 
stant Chf=i?hf /m ~ 15T/^s, and in a number of other 
intermetallics with similar Chf value. In the 122 family of 
isostructural pnictides, the ratios of hyperfine field and 
moment are smaller and present rather scattered values 
(see Table HH where the moment is taken from neutron 
diffraction data). A strict proportionality between mo- 
ment and hyperfine field can be however questioned in the 
pnictides. Indeed, the magnetic susceptibility has been 
found to be anisotropic (xx/X// — 1-5-2) both in un- 
doped BaFe2As2 [H[ and in Co-doped compounds [23| . 
This is a hint to the presence of an "unquenched" or- 
bital moment at the Fe site. This implies a probably 
small, but non-zero contribution of the d-orbitals to the 
hyperfine field, with a sign opposite to that of the Fermi 
contact field. Then, like in Fe 2+ compounds (24[, the 
hyperfine field is no longer strictly proportional to the 
moment. The hyperfine field modulation then would not 
reflect exactly that of the Fe moment. 

However, assuming as a first good approximation that 
proportionality holds in the Ba(Fei_ a; Co :c )2As2 series 
with the same constant as in BaFe2As2, then the max- 
imum Fe moment of the modulation at 4.2 K would be 
0.96, 0.84 and 0.72 fi B for 1=0.014, 0.03 and £=0.045 
respectively. 



D. Thermal variation of the moment and nature of 
the SDW transition 

In the neutron diffraction studies of the series 
Ba(Fei_ K Co2;)2As2 [a. l9Hllj. the integrated intensity of 
the magnetic Bragg peaks continuously drops towards 
zero as temperature is increased towards Tsdw- This 
behaviour is in contrast with the 75 As NMR work in 
BaFe2As2 0, where the transferred hyperfine field at 
the As nucleus shows a weak thermal variation and it 
was concluded that the SDW transition is first order. A 
similar discrepancy is observed in the present work in 
Co-doped materials, where the main Fourier component 
hi of the modulation retains a sizeable value at Tsdw 
(see FigJB]Top). However, a correct comparison with the 
neutron data must take into account the fact that the 
intensity of a magnetic Bragg peak is proportional to the 
square of the magnetic moment. In case of a modulated 
structure with Fourier coefficients m2fe+i, it should scale 
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FIG. 8: [color online] Variation of the normalised values of the 
Bragg peak intensity for x=0.04 from Ref. [l^l (diamonds) and 
of the Mossbauer derived {m 2 (x,T)) for 1=0.03 from Eqnf2] 
(squares), as a function of the reduced temperature T/Tsdw- 
The Mossbauer data above Tsdw belong to the short range 
AF regions. 



with (m 2 ), where the brackets denote an average over a 
period of the moment modulation: 



(m 2 (x,T)) 



E 5 

k=0,n 



'2k- 



-iO,T) cx i h ik+i(x,T). 



k=Q,n 



(2) 

Figure[8]shows the thermal variations of (m 2 (x, T)) (from 
the Ji2k+i(x,T) values determined by Mossbauer spec- 
troscopy) for £=0.03 and of the intensity of a mag- 
netic Bragg peak for £=0.04 [l^. The (m 2 ) value is 
reduced by half at Tsdw, m disagreement with the neu- 
tron data. Similar behaviours occur for x=0.014 and 
0.045. It can thus be concluded that the SDW transi- 
tion in the Ba(Fei_2; Co^ As2 series, as observed with 
the Mossbauer local probe, also presents a strong first 
order character. 

Comparing the moment values obtained from the neu- 
tron data [llf for x=0.025 and 0.045 with the Mossbauer 
results is more difficult since, in the scaling of the neutron 
Bragg intensity, it was assumed that the magnetic struc- 
ture is identical to that of BaFe2As2, i.e. commensurate 
with the lattice spacings. The mean saturated Fe mo- 
ment value m s obtained from neutron diffraction should 
be compared with: 



m s (x) = v/(m 2 (x,T = 4.2 K)). 



(3) 



Assuming proportionality between hyperfine field and 
moment, and using the Mossbauer derived values for 
h 2k +i(x,T = 4.2 K), one gets: m s =0.70, 0.59 and 
0.45 hb resp. for a;=0.014, 0.03 and 0.045. This must be 
compared with m s =0.35 and 0.17/zb resp. for £=0.025 
and 0.045 derived from the neutron data in Ref . (TTI| . 

Further work is needed to resolve the disagreement be- 
tween local probe and neutron diffraction data about the 
Fe moment variation with temperature and Co doping. 
VI. CONCLUSIONS 

57 Fe Mossbauer absorption spectra have been mea- 
sured in the Ba(Fei_ x Co x )2As2 series for x=0.014, 0.03 
and 0.045. In all the samples, unusual magnetic hyper- 
fine lineshapes are observed from 4.2 K to Tsdw- They 
are interpreted in terms of a continuous modulation of 
hyperfine fields at the Fe sites arising from an incom- 
mensurate electronic SDW state. Incommensurability is 
present at the lowest doping level, i.e. for £=0.014, while 
the magnetic structure in BaFe2As2 is commensurate. 
Above Tsdw j we do not observe a purely paramagnetic 
signal, but a coexistence of magnetic hyperfine and para- 
magnetic subspectra in a large temperature range, up to 
Taf ~ 1-5 Tsdw approximately for all samples. We be- 
lieve this is due to the persistence of AF short range cor- 
related regions, probably dynamic, well above the SDW 
transition, which is found to have a strong first order 
character. Our analysis of the hyperfine field modula- 
tions as a function of x shows that magnetism is not 
strongly depleted as doping increases, i.e. the mean Fe 
moment or the maximum of the modulation does not de- 
crease with x as rapidly as inferred from neutron diffrac- 
tion data. We believe the incommensurability of the mag- 
netic structure has not been evidenced by neutron scat- 
tering up to now due to the probable very small value 
of the incommensurate component (s) of the SDW wave 
vector. This makes Mossbauer spectroscopy a valuable 
tool to investigate these structures: although the prop- 
agation vector cannot be obtained with this technique, 
the spectral shape is sensitive to minute incommensura- 
bilities. 
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